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Abstract 

Measurements are presented of diffractive open charm production at HERA. The event 
topology is given by ep — ► eXY where the system X contains at least one charmed hadron 
and is well separated by a large rapidity gap from a leading low-mass proton remnant 
system Y. Two analysis techniques are used for the cross section measurements. In the first, 
the charm quark is tagged by the reconstruction of a Z)*^(2010) meson. This technique is 
used in deep-inelastic scattering (DIS) and photoproduction (jp). In the second, a method 
based on the displacement of tracks from the primary vertex is used to measure the open 
charm contribution to the inclusive diffractive cross section in DIS. The measurements are 
compared with next-to-leading order QCD predictions based on diffractive parton density 
functions previously obtained from a QCD analysis of the inclusive diffractive cross section 
at HI. A good agreement is observed in the full kinematic regime, which supports the 
validity of QCD factorization for open charm production in diffractive DIS and *yp. 
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1 Introduction 



Diffractive processes in positron-proton (ep) collisions are those where the hadronic final state 
is separated by a large gap in rapidity, without hadrons, into two systems X and Y, where the 
system Y may consist only of a proton or low mass system. The system X is known as the 
photon dissociative system. The diffractive event signature is understood to arise from a color 
singlet exchange between the two systems X and Y. 

In quantum chromodynamics (QCD), the theory of strong interactions, the hard scattering 
collinear factorization theorem [1] predicts that the cross section for diffractive deep-inelastic ep 
scattering (DIS) factorizes into a set of universal diffractive parton density functions (DPDFs) 
of the proton and process-dependent hard scattering coefficients. Next-to-leading order (NLO) 
DPDFs have been determined by QCD fits to the measured cross sections of inclusive diffractive 
scattering at HERA [2, 3] within the factorizable Pomeron model [4] and using the DGLAP 
evolution equations [5]. The DPDFs have been found to be dominated by gluons, which carry 
~ 70 % of the momentum of the diffractive exchange. 

If QCD factorization is fulfilled, NLO QCD calculations based on the diffractive parton 
density functions of [2, 3], should be able to predict the production rates of more exclusive 
diffractive processes like dijet and open charm production in shape and normalization. For 
diffractive dijet production this has been tested in photoproduction {^p) and in DIS [6]. In 
the regime of DIS the predictions of QCD have been found to be in good agreement with the 
experimental results. 




Figure 1: The main processes of diffractive open charm production at HERA in the collinear 
factorization approach. Figure a) shows the direct process where the photon enters the hard 
scatter itself. Figure b) shows the resolved photon process where only a reduced fraction x 7 of 
the photon 's momentum takes part in the hard scatter. 



In the collinear factorization approach diffractive open charm production at HERA is ex- 
pected mainly to proceed via boson gluon fusion (BGF) as depicted in Figure [T^. Thus it is 
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directly sensitive to the gluon content of the diffractive exchange, which is only determined 
indirectly and for low momentum fractions zjp of the gluon in inclusive diffractive scattering 
via scaling violations [3]. In the BGF process a charm quark anti-quark pair (cc) is produced of 
which one quark couples to the photon with virtuality Q 2 and the other to a gluon that emerges 
from the diffractive exchange. 

In Figure [T^ the "direct photon" process is shown, where the photon itself enters the hard 
scatter, which is expected to be dominant for cc-production in DIS and photoproduction. In 
photoproduction, however, the quasi real photon may also evolve into a hadronic structure, as 
indicated in Figure before it enters the hard scatter. In this case only a fraction x 7 < 1 
of the photon's momentum takes part in the scattering process; the rest forming a remnant. In 
these "resolved photon" processes initial state interactions may take place between the photon- 
and the proton-remnant systems, destroying the rapidity gap signature and thus the diffractive 
nature of the process. A breakdown of QCD factorization has been observed for diffractive dijet 
production in pp collisions at the Tevatron [7], where the prediction overestimates the observed 
rate by approximately one order of magnitude. Diffractive open charm production is especially 
suitable for testing a potential suppression of the direct photon component of the production 
mechanism in photoproduction. 

In an alternative theoretical approach DPDFs are not introduced and diffractive scattering is 
explicitly modeled by the perturbative exchange of a colorless gluon state (two gluons or a gluon 
ladder). Formulated in the proton rest frame the "two-gluon" state of the proton can couple 
directly to the cc pair (YP ~~ * °cp) or to a ccg color dipole fluctuation of the photon (7*^ — > 
ccgp) [8]. The gluon density of the proton is usually determined from fits to the inclusive DIS 
cross section in the A; r factorization [9] scheme. A model combining the perturbative two-gluon 
approach with the collinear factorization scheme, which has also been used to fit the HERA 
diffractive DIS cross sections, is given in [10]. 

Two methods to identify charm production are presented in this paper. In the first method 
the charm quark is tagged by the reconstruction of D* mesons. The measurement is performed 
in DIS and, due to the high selectivity of the trigger, extended to jp, where it represents the first 
cross section measurement of diffractive open charm production at HERA. In DIS it supercedes 
a former analysis of HI [11] with increased statistics and with reduced systematic uncertainties. 
A similar measurement in DIS was performed by the ZEUS collaboration [12]. The results are 
presented in the form of integrated and differential D* cross sections and in DIS are extrapolated 
into the unmeasured phase space of the D* meson using NLO QCD calculations in order to 
determine the open charm contribution to the diffractive cross section. The second method, 
which was used to measure the total inclusive charm and beauty cross sections in DIS [13, 14], 
is used here for the first time in diffractive DIS. In this method, referred to in the following as the 
'displaced track analysis', the charm quark is identified by the reconstruction of tracks, which 
are displaced from the interaction vertex, that arise due to long lived charmed hadrons. This 
method is used in a kinematic region with high acceptance for the decay products of charmed 
hadrons within the silicon vertex detector of HI, which is used in the reconstruction of these 
tracks. With this method it is thus possible to measure the total open charm contribution to the 
diffractive cross section with small extrapolations from QCD calculations. 

In section |2] the kinematic variables used throughout the paper are introduced. A short 
discussion of the HI detector and the event selection are given in sections |3] and |4j followed by 



5 



a description of the event simulation in section |5] and the cross section determination with the 
two independent methods in section The NLO QCD calculations and the comparison of the 
measured cross sections with the NLO QCD calculations are discussed in sections |7] and |5] 



2 Kinematics of Diffractive ep Scattering 

Due to the diffractive nature of the process the photon (with four-momentum q) and the proton 
(with four-momentum P) dissociate into two distinct hadronic systems X and Y (with four- 
momenta px and py, respectively), which are separated by a large gap in rapidity between the 
final state hadrons. The kinematics of the inclusive ep scattering are fully determined by the 
negative squared four momentum transfer of the exchanged photon Q 2 , the squared center of 
mass energy of the ep scattering process s and the inelasticity y. In addition, the following 
variables are defined to characterize the diffractive nature of the process 



M Y =p Y ; t=(P-p Y ) ; X]p = ; P = - — r, (1) 

q-P 2q-(P-p Y ) 

where My and t denote the invariant mass of the system Y and the squared four-momentum 
transferred at the proton vertex, respectively. The variable x p can be interpreted as the longitu- 
dinal momentum fraction of the diffractive exchange with respect to the proton. The variable j3 
(which is only defined in DIS) corresponds to the Bjorken x variable from inclusive scattering 
taken with respect to the diffractive exchange. The quantities t and M Y are constrained to be 
small by the experimental selection and are integrated over implicitly. For the D* analyses the 
observable Zjp is introduced as 

4- - ^ ( 2) 

xp-y ■ s 

where s obs is a hadron level estimate of the invariant mass of the cc pair emerging from the 
hard scattering process. It is reconstructed from the scattered positron and the kinematics of 
the reconstructed D* meson including an approximate correction of the momentum of the D* 
meson to the momentum of the charm quark [11]. In direct BGF processes z°p is a direct 
estimator for the longitudinal momentum fraction zjp of the gluon that enters the scattering 
process with respect to the momentum of the diffractive exchange. In resolved processes zp 
cannot be disentangled by the reconstruction method from the momentum fraction x 7 that enters 
the hard scattering process from the photon side. 



3 The HI Detector 

A detailed description of the HI detector can be found in [15]. Only the components most 
relevant for this analysis are briefly discussed here. The coordinate system is centered at the 
nominal ep interaction point with the z-axis pointing along the beam direction of the outgoing 
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proton, also referred to as the 'forward' direction in the following. Charged particles emerging 
from the interaction region are measured by the Central Tracking Detector (CTD), which covers 
a range of —1.74 < 77 < 1.74 in pseudorapidity l . The CTD comprises two large cylindrical 
Central Jet drift Chambers (CJCs) and two z chambers situated concentrically around the beam- 
line within a solenoidal magnetic field of 1.15 T. It also provides triggering information based 
on track segments measured in the r-0 plane of the CJCs, and on the z position of the event 
vertex obtained from the double layers of two multi-wire proportional chambers (MWPCs). The 
CTD tracks are linked to hits in the central silicon tracker (CST) [16] to provide precise spatial 
track reconstruction. The CST consists of two layers of double-sided silicon strip detectors 
surrounding the beam pipe, with a coverage of —1.3 < rj < 1.3 in pseudorapidity for tracks 
passing through both layers. 

The tracking detectors are surrounded by a Liquid Argon calorimeter (LAr) in the forward 
and central region (—1.5 < rj < 3.4) and by a lead- scintillating fiber calorimeter (SpaCal) 
with electromagnetic and hadronic sections in the backward region [17] (—4 < i] < —1.4). 
These calorimeters provide energy and angular reconstruction for final state particles from the 
hadronic system. DIS events are identified by the energy deposits of the scattered positron in 
the SpaCal calorimeter. Photoproduction events are selected with a crystal Cerenkov calorime- 
ter located close to the beam pipe at z = —33.4 m in the positron direction (electron tagger), 
which measures the energy deposits of positrons scattered by angles of less than 5 mrad. An- 
other Cerenkov calorimeter located at z = — 103 m is used to determine the ep luminosity by 
detecting the radiated photon emitted in the Bethe-Heitler process (ep — > epj). 

For the rapidity gap selection a set of detectors close to the beam pipe in the forward direc- 
tion is used. The Forward Muon Detector (FMD) is located at z = 6.5 m and covers a pseudora- 
pidity range of 1.9 < rj < 3.7. It may also detect particles produced at larger 77 due to secondary 
scattering within the beam pipe. A PLUG hadronic sampling calorimeter allows energy mea- 
surements in the range of 3.5 < rj < 5.5. Finally, particles in the region of 6.0 < rj < 7.5 can 
be detected by the Proton Remnant Tagger (PRT), a set of scintillation counters surrounding the 
beam pipe at z = 26 m. 



4 Event Selection 

The data presented in this analysis were collected over the years 1999 and 2000 and corre- 
spond to an integrated luminosity of C = 47.0 pb _1 for the D* analyses and 48.3 pb" 1 for the 
displaced track analysis. At this time HERA was operated with positrons of energy 27.6 GeV 
and protons of energy 920 GeV so that the center of mass energy of the ep collision is y/s = 
318 GeV. 

DIS events are triggered by an electromagnetic energy cluster in the SpaCal calorimeter. 
In the D* analyses the trigger further requires a charged track signal in the CTD and a recon- 
structed event vertex, while a looser track requirement of hits in the MWPCs is used in the 
trigger for the displaced track analysis. In the offline analyses the scattered positron is selected 
as an electromagnetic SpaCal cluster with energy E' e > 8 GeV. Photoproduction events are 

'The pseudorapidity r\ of an object detected with polar angle 6 is defined as rj = — In tan(0/2). 
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suppressed by requiring ~ Pz,%) > 35 GeV. Here, E t and p Zji denote the energy and 

longitudinal momentum components of a particle and the sum is over all final state particles 
including the scattered positron and the hadronic final state (HFS). The HFS particles are re- 
constructed using a combination of tracks and calorimeter deposits in an energy flow algorithm 
that avoids double counting. The z position of the interaction vertex is required to lie within 
±35 cm (±20 cm) of the center of the CTD for the D* (displaced track) analyses, where the 
reduced range in the displaced track analysis is chosen in order to match the smaller acceptance 
of the CST The kinematic variables of the DIS scattering process Q 2 and y are reconstructed 
using a method which uses the angle of the hadronic final state in addition to the energy and the 
polar angle of the scattered positron [3]. The accepted kinematic range in DIS is restricted to 
2 < Q 2 < 100 GeV 2 and 0.05 < y < 0.7 for the D* analysis and to 15 < Q 2 < 100 GeV 2 
and 0.07 < y < 0.7 for the displaced track analysis, where the reduced kinematic range in 
the displaced track analysis is chosen such that the direction of the quark struck by the pho- 
ton mostly lies within the angular acceptance of the CST and that the HFS has a significant 
transverse momentum. 

Photoproduction events are selected by a trigger that requires a scattered positron to be 
measured in the electron tagger, a charged track signal in the CTD and a reconstructed event 
vertex. The events have passed an additional online software filter that selects events with 
candidates for charmed hadron decays by calculating the invariant mass of track combinations. 
The inelasticity y is reconstructed from the energy of the scattered positron and is restricted to 
the range 0.3 < y < 0.65. The photon virtuality is experimentally restricted to Q 2 < 0.01 GeV 2 . 

In all analyses presented in this paper diffractive events are selected by the absence of 
hadronic activity above noise thresholds in the most forward part of the LAr calorimeter (77 > 
3.2) and in the forward detectors. This selection ensures that the gap between the systems X 
and Y spans more than 4 units between rj = 3.2 and 7.5 in pseudorapidity. As My is not di- 
rectly measurable by this method the data are corrected to a visible range of My < 1.6 GeV 
and \t\ < 1 GeV 2 , consistent with former measurements [3,6, 1 1], with the help of Monte Carlo 
simulations. The variable x F is calculated from 



= M 2 X = J2 (El - Pi - p\ y - Pi) , (3) 

s-y i 

where the sum for the calculation of Mx runs over all HFS objects in the system X. Each of 
the presented analyses is restricted to xp < 0.04, which suppresses contributions from non- 
diffractive scattering and secondary Reggeon exchanges. The displaced track analysis is further 
restricted to M x > 6 GeV. 



5 Event Simulation and Acceptance Correction 

The data are corrected for trigger efficiencies, detector acceptances, efficiencies, and migration 
effects due to the finite resolution of the HI detector using a Monte Carlo simulation. All the 
generated events are passed through a detailed simulation of the detector response based on the 
GEANT simulation program [18] and reconstructed using the same reconstruction software as 
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used for the data. For the event simulation residual noise contributions in the LAr calorimeter 
and the forward detectors are taken into account. 

Events are generated using the RAPGAP event generator [19], which simulates the process 
e + p — > e + Xp with xp < 0.15, assuming proton vertex factorization. Both Pomeron and 
Reggeon sub-leading exchanges are included. The t dependence is of the form da/dt oc e BEht 
with a slope parameter B EL = 6 GeV -2 . For the simulation of diffractive events containing 
charm quarks RAPGAP implements the BGF process in leading order (LO) of pQCD. For the 
D* analyses LO DPDFs are taken from a former analysis of HI [2]. For the displaced track 
analysis the DPDFs are taken from [3]. To simulate higher order effects of QCD, parton show- 
ers are included in the calculations. Fragmentation is performed according to the Lund string 
model [20]. In DIS RAPGAP is interfaced to the QED simulation program HERACLES [21] to 
evaluate the radiative effects of QED. For diffractive photoproduction the contributing diagrams 
of charm excitation and other resolved photon processes are included in the event generation, 
using the LO parton distribution functions for the resolved photon obtained in [22]. In the exci- 
tation processes the charm quark is treated as a massless parton in the resolved photon, whereas 
in all other processes the charm mass is taken into account in the calculations. The resolved pro- 
cesses are found to contribute less then 10% of the charm signal and to be mainly concentrated 
at large values of xjp and small values of pt(D*) and Zp s . 

Due to the limited detector acceptance in the forward region of HI it is not possible to 
efficiently detect a break-up of the proton into a low mass resonant state Y. To keep the un- 
certainties arising from such proton dissociation processes small the measurement is integrated 
over the region M Y < 1.6 GeV and \t\ < 1 GeV 2 . Diffractive proton dissociative events 
in the region M Y < 5 GeV are simulated using RAPGAP with a cross section dependence 
of the form e BpDt with B PD = 1.6 GeV~ 2 and an approximate M Y dependence of the form 
da/dMy oc 1/M Y [23]. The correction factor S pdls for migrations across the measurement 
boundary is evaluated in the simulation for each kinematic bin. In the simulations the ratio of 
proton elastic to proton dissociative interactions is taken to be 1 : 1, which is in accordance with 
the inclusive measurements of [3, 24]. The value of 5 pdm is found to be in the range 0.88 — 0.97. 

Non-diffractive events with M Y > 5 GeV or xp > 0.15 are simulated by RAPGAP in DIS 
and by the event generator pythia [25] in photoproduction. The non-diffractive background 
contribution in the final event selections is estimated to be less than 3% for all data samples. 



6 Open Charm Selection 

Charm quarks are selected by two independent methods. In the first method they are selected 
by the full reconstruction of D* mesons. This provides a clear signature, which enables the 
tagging of charm quarks in DIS and photoproduction. In the second method the more general 
character of the long lifetime of charmed hadrons is used, by reconstructing the displacement 
of tracks from the primary vertex in the CST of HI, similarly to inclusive charm production 
measurements in [13, 14]. This provides the advantage of a high acceptance for charm quarks 
and small correction factors for extrapolations to the full phase space. It is therefore especially 
suited for a measurement of the total diffractive charm cross section. 
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6.1 Diffractive D* Analyses 



In the D* analyses D*^ 1 mesons are fully reconstructed using the decay channel 

D* + - D\+ ow -> (K-n+)n+ ow (+C.C), (4) 

which has a branching ratio of 2.57% [26]. The decay products are detected in the CTD. To 
ensure good detection efficiency and to reduce combinatorial background, the tracks are re- 
quired to lie within an angular range of 20° < 9 < 160° and to have a transverse momentum 
p t relative to the beam axis of at least 120 MeV for the ix s i ow , 300 MeV for the other pion 
and 500 MeV for the K candidate. The invariant mass of the Kix combination has to be con- 
sistent with the D° mass within ±80 MeV. The transverse momentum and pseudorapidity of 
the reconstructed D* meson candidate are restricted to pt{D*) > 2 GeV and |^(D*)| < 1.5. 
The distributions of the mass difference AM = M{K^TT ±r i:f low ) — M^K^n^) for all track 
combinations which fulfill the above requirements for all selected events in DIS and photo- 
production are shown in Figure El The number of D*^ mesons is determined by fitting these 
distributions with a Gaussian function for the signal plus a background parameterization given 
by N(AM — m 7r ) Ue (l — u s (AM) 2 ), where m w denotes the mass of the charged pion and N, u e 
and u s are free parameters. The position and the width of the Gaussian function are fixed to val- 
ues taken from higher statistics samples where no diffractive cuts were applied. The resulting 
numbers of identified D*^ 1 mesons in DIS and photoproduction are summarized in table [lj 

Differential cross sections are obtained from the fitted number of D* mesons in each mea- 
surement bin. A correction is applied for mass reflections originating from decays of the D° 
meson other than that given in equation|4j which has been estimated to be 3.5% of the D* sig- 
nal [27]. A correction factor of ~ 0.95 for the effects of initial and final state QED radiation 
is applied to the DIS cross sections. The cross sections are bin center corrected using RAPGAP 
to determine the point in the bin at which the bin-averaged cross section equals the differential 
cross section. 

6.2 Displaced Track Analysis 

The production of open charm in diffraction is also investigated using a largely independent 
method, which has been used in [13] and [14] to measure the total inclusive charm and beauty 
cross sections in DIS. This method distinguishes events containing heavy quarks from those 
containing only light quarks by reconstructing the displacement of tracks from the primary 
vertex in the transverse plane (impact parameter), caused by the long lifetimes of the charm and 
beauty flavored hadrons, using the precise spatial information from the CST of HI. Due to the 
low beauty fraction in the diffractive data sample, it is not possible to make a measurement of 
the beauty cross section and only a measurement of the charm cross section is presented in this 
paper. 

As in [13, 14] the primary event vertex in the r-0 plane is reconstructed from all tracks 
(with or without hits in the CST) using the information on the position and transverse extent 
of the beam interaction region. For the analysis, tracks are selected if they have a transverse 
momentum of more than 0.5 GeV and at least two associated hits in the CST. The impact 
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parameter of a track is defined as the distance of closest approach (DCA) of the track to the 
primary vertex point in the transverse plane. 

In order to determine a signed impact parameter (5) for a track, the azimuthal angle of 
the struck quark quar k must be determined for each event. To do this, jets with a minimum 
p t of 2.5 GeV, in the angular range 15° < 9 < 155°, are reconstructed using the invariant 
k t algorithm [28] in the laboratory frame using all reconstructed HFS particles. The angle 
0quark is defined as the of the jet with the highest transverse momentum or, if there is no jet 
reconstructed in the event, as 180° — e iec> where clcc is the azimuthal angle of the scattered 
positron in degrees. Monte Carlo simulations indicate that ~ 78% of all charm events have at 
least one reconstructed jet in the kinematic region described above. The direction defined by the 
primary vertex and qua rk in the transverse plane is called the 'quark axis'. If the angle between 
the quark axis and the line joining the primary vertex to the point of DCA of the track is less 
than 90°, 5 is defined as positive, and is defined as negative otherwise. Tracks with azimuthal 
angle outside ±90° of qU ark are rejected. The estimated error on § is denoted as a(S). 

To distinguish between the charm and light quark flavors a similar method to that in [14] is 
used. The quantity Si (S2) is defined as the significance (8/a(5)) of the track with the highest 
(second highest) absolute significance that is associated to the quark axis. In the present analysis 
S3, which is the significance of the track with the third highest absolute significance, is not used 
due to lower statistics than in [14]. Events where Si and S2 have opposite signs are excluded 
from the S2 distribution, but contribute to the S\ distribution. The distributions of Si and 5*2 are 
shown in Figure |3] for the kinematic region given in section UJ A reasonable description of the 
data by the simulation is observed. The light quark significance distributions are approximately 
symmetric around zero, whereas the charm distributions have an excess in the positive bins 
compared with the negative. It is thus possible to substantially reduce the uncertainty due to 
the resolution of 5 and the light quark normalization, by subtracting the contents of the negative 
bins in the significance distributions from the contents of the corresponding positive bins. The 
subtracted distributions are shown in Figure |U 

The fractions of charm and light quark flavors in the data are extracted in three M\ inter- 
vals using a least squares simultaneous fit to the subtracted S\ and Si distributions (as shown 
in Figure |4j) and the total number of reconstructed diffractive events before any track selection. 
The significance distributions of the charm, beauty and light flavors, as predicted by the Monte 
Carlo simulation for the luminosity of the data, are used as templates. In each interval the charm 
and light quark flavor contributions from the Monte Carlo simulation are scaled by factors P c 
and Pi, respectively, to give the best fit to the observed subtracted S\, S2 distributions and the 
total number of events. Since the same event may enter the Si and the S2 distributions, it was 
checked using a high statistics Monte Carlo simulation that this has a negligible effect on the 
results of the fits with the statistics of the present data. Only the statistical errors of the data and 
the Monte Carlo simulation are taken into account in the fit. The beauty scale factor is fixed to 
Pb = 1, and varied in the evaluation of the systematic uncertainties (see section IQl) . The results 
of the fit to the complete data sample are shown in Figure |4j The fit gives a good description 
of all significance distributions, with a y 2 /n.d.f of 18.1/12. Values of P c = 0.77 ± 0.09 and 
Pi = 0.97 ± 0.03 are obtained. It can be seen that the resulting distributions are dominated by 
charm quark events, the light quarks contributing only a small fraction, mainly due to strange 
hadrons, for all values of the significance. The beauty contribution forms a small fraction over- 
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all, but increases with increasing M x . Acceptable x 2 values are also found for the fits to the 
samples in the separate M\ intervals. 

The results of the fit in each M x interval are converted to a measurement of the diffractive 
differential cross section using: 

dVf _ p c N™ c z™5™ d 5 pdis 5f cc 
dQ 2 dydM 2 x ~ £-BV ' ( } 

where A^ ( ! vlCgon is the number of generated charm events expected from the Monte Carlo simu- 
lation in each bin with volume BV corresponding to the luminosity of the data C. A bin center 
correction 5® cc in the range 0.89 — 1.21 is calculated using the NLO QCD expectation to cor- 
rect the bin averaged cross section to the cross section at a specified point in Q 2 , y and M\. 
A correction factor of S r c ad ~ 0.93 for initial and final state QED radiation is applied. The 
correction factor for proton dissociation S pdls is described in section|3J 

Measurements of the ratio of the diffractive charm cross section to the total diffractive cross 
section are made where the total diffractive cross section is determined using 

d 3 a D jyrecjyMCgen^rad^pdis^BCC 

(6) 



dQ 2 dy dM 2 x n MCicc C ■ BV 

where iV rcc is the number of reconstructed data events in the bin after the event selection de- 
scribed in section HI jV MCrec (jV MCgen ) is the number of reconstructed (generated) Monte Carlo 
events in the bin; 5 rad and <5 BCC are the radiative correction and bin center correction for inclu- 
sive diffraction, respectively. The ratio is then given by: 



f cc = dVg d 3 a D 
Jd dQ 2 dy dM 2 x 1 dQ 2 dy dM 2 x ' 



6.3 Systematic Uncertainties 

The following sources of systematic uncertainty for the two different analysis methods are taken 
into account; the estimated values are given in tabled 

• The simulated trigger efficiencies for the D* analyses are compared with the efficiencies 
determined from data using monitor trigger samples. Within the statistics of these data 
samples the simulated trigger efficiencies are found to agree with the data, with a remain- 
ing uncertainty in the range ±(3 — 5)% depending on the analysis. For the displaced track 
analysis an uncertainty of 1% is assigned as determined from the data. 

• For the DIS measurements the reconstructed polar angle and the energy of the scattered 
positron are varied within the estimated uncertainties of ±1 mrad for the angular mea- 
surement and ±1% for the energy scale of the SpaCal, leading to an uncertainty of ±2% 
on the cross section measurements. In photoproduction a variation within the estimated 
uncertainty of ±1.5% on the energy scale of the crystal Cerenkov calorimeter of the elec- 
tron tagger results in an uncertainty of ±2%. 
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The uncertainty of the track reconstruction efficiency and uncertainties related to the sig- 
nal extraction for the D* analyses have been determined by analyzing inclusive D* sam- 
ples as in [29] and are estimated to be ±6% for the reconstruction efficiency of the three 
daughter tracks of the D* meson in the CTD and ±6% for the signal extraction. The 
uncertainty on the correction for mass reflections is estimated to be ±1.5% [27]. 

For the displaced track analysis a track efficiency uncertainty of ±2% due to the CTD 
and of ±1% due to the CST is estimated, resulting in an uncertainty of ±2% on the 
cross sections. An uncertainty in the resolution of 5 of the tracks is estimated by varying 
the resolution by an amount that encompasses the differences between the data and the 
simulation (see Figure |3]). This is achieved by applying an additional Gaussian smearing 
in the Monte Carlo simulation of ±200 /im to 5% of randomly selected tracks and ±25/im 
to the rest, resulting in an error of 2% on the cross sections. 

The effect of a ±4% uncertainty in the energy scale of the hadronic final state leads to a 
change of the cross section in the range ±(1 — 3)% depending on the analysis. 

The uncertainty in the acceptance and migration corrections due to uncertainties in the 
physics models for diffractive charm production is estimated by varying the shape of var- 
ious kinematic distributions in the Monte Carlo simulation within limits set by the present 
measurements. Reweighting the shapes of the xjp, (3 and Q 2 distributions by (:^) ±0 ' 25 > 
^±0.3 and ^ _|_ i ogl0 (g 2 /GeV 2 )) ±1 in DIS results in an uncertainty of ±5% on the total 
cross section for the D* analysis and ±(12 — 18)% for the displaced track analysis. A 
variation of the xjp and y distributions by (^■) ±0 ' 25 and (^) ±0 ' 2 in photoproduction re- 
sults in an uncertainty of ±1% on the total cross section. The uncertainty on the fraction 
of the Reggeon contribution is estimated by varying its normalization in the simulation 
by ±100%, which leads to an uncertainty of ±1% (±4%) for the D* analyses in DIS 
(photoproduction) and ±(1 — 9)% for the measurement bins of the displaced track analy- 
sis. A variation of the t distribution by e ±2t for proton elastic scattering and the M Y and 
the t distribution by (^r) ±0 ' 3 and e ±11 for proton dissociative scattering as well as the 
ratio of proton elastic to proton dissociative scattering between 1 : 2 and 2 : 1 results in a 
systematic uncertainty on the cross sections in the range ±(4 — 5)%. 

The uncertainties for residual noise in the FMD and the PLUG calorimeter in the simula- 
tion are estimated using a set of randomly triggered events during the data taking period 
and result in a combined uncertainty of ±1.5%. The tagging efficiency of the PRT for 
proton dissociative systems with My > 1.6 GeV or \t\ > 1 GeV 2 in the simulation is 
adjusted with the help of an independent non-diffractive data sample with activity in the 
forward part of the LAr calorimeter and the FMD, where such events are enriched. The 
effect of the remaining uncertainty on this efficiency on the cross section measurements 
is estimated by varying the simulation within the statistical accuracy of the measured 
efficiency and is estimated to lie between ±(7 — 9)%. The uncertainty on the tagging effi- 
ciency of the FMD is estimated to be ±10% [3]. The effect on the cross sections is ±1%. 
The residual influence of non-diffractive background from events without a rapidity gap 
is estimated by assigning a ±100% uncertainty to the corresponding event samples in 
the RAPGAP simulation. This leads to an uncertainty on the cross sections in the range 
±(1 - 3)%. 
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• The uncertainty of the charm fragmentation scheme is estimated by changing the para- 
metrization of the longitudinal fragmentation function from the Lund-Bowler model [20] 
to Peterson functions with e = 0.078 (e = 0.058) [30] in the simulation of the events for 
the D* (displaced track) analyses, which results in an uncertainty on the cross section of 
±1% (±4%) for the D* analyses in DIS (photoproduction) and of ±7% in the displaced 
track analysis. 

• For the displaced track analysis the uncertainties on the lifetimes of the various D mesons, 
decay branching fractions and mean charge multiplicities are estimated by varying the 
input values of the Monte Carlo simulation by the errors on the world average measure- 
ments. The values and the uncertainties for the lifetimes of the D mesons are taken 
from [26] and those from the branching fractions of charm quarks to hadrons from [31]. 
They are consistent with measurements in DIS at HERA [32]. The values and the uncer- 
tainties for the mean charged track multiplicities for charm quarks are taken from [33]. A 
combination of all these uncertainties results in an error of 3% on the cross sections. For 
the D* analyses the uncertainty of ±2.5% on the branching fraction for the decay channel 
in equation|4]is taken from [26]. 

• The uncertainty on the asymmetry of the 5 distribution for the light quarks in the dis- 
placed track analysis is estimated by repeating the fits with the subtracted light quark 
significance distributions (shown in Figure |U) changed by ±50%. The light quark asym- 
metry is checked to be within this uncertainty by comparing the asymmetry of the Monte 
Carlo events to that of the data, in the region of 0.1 < \5\ < 0.5 cm, where the light quark 
asymmetry is enhanced. This results in an uncertainty on the cross section of ±4% at 
high M x and of ±16% at low M x . 

The uncertainty on the beauty contribution for the displaced track analysis is estimated 
by repeating the fits with the subtracted beauty quark significance distributions (shown in 
Figure |4j) changed by tioo%, which results in an negligible error on the cross section at 
low Mx increasing to ^5 % and +13% in the middle and high Mx bins, respectively. 

• An uncertainty on the quark axis in the displaced track analysis is estimated by shifting 
it by ±2° (±5°) for events with (without) a reconstructed jet. These shifts have been 
estimated in [14] by comparing the difference between quar k and the track azimuthal 
angle in data and Monte Carlo simulation. The resulting error on the cross sections is 
±3%. 

• The uncertainty in the calculation of QED radiative effects is found to be ±2% in DIS. 

• The uncertainty in the bin center correction for the displaced track analysis is estimated 
by varying the shape of the Q 2 , (3 and xjp distributions of the NLO QCD expectation. 
This leads to a ±(8 — 10)% uncertainty on the cross sections. 

• The uncertainty of the luminosity determination is estimated to be ±1.5%. 

The total systematic uncertainty for each data point has been obtained by adding all indi- 
vidual contributions in quadrature. For the D* analyses it ranges between 15% and 30% for the 
differential cross sections and amounts to ±15% for the integrated cross section in both kine- 
matic regimes. For the displaced track analysis it ranges between 26% and 47% for the three 
points of the inclusive charm cross section measurement. 
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7 QCD Calculations 



7.1 NLO Calculations in Collinear Factorization 

The measured charm cross sections are compared with NLO QCD calculations based on two 
alternative sets of diffractive parton density functions from HI [3] which both provide a good 
description of the inclusive diffractive DIS data. As default the standard parameterization HI 
2006 DPDF Fit A is chosen. The alternative set of DPDFs (HI 2006 DPDF Fit B) is obtained 
from a slightly different parameterization of the gluon density at the starting scale of the fit 
procedure. It leads to a steeper fall-off of the gluon density at higher values of z F . In the fit to 
the inclusive diffractive DIS data [3] charm quarks are treated as massive, appearing via BGF- 
type processes up to order a 2 s [34]. The quark mass is set to m c — 1.4 GeV and the scale for 
heavy flavor production to /i r — fif — 2m c . 

In order to be able to compare the measured D* cross section to the results based on the 
NLO QCD fits diffractive versions of the programs HVQDIS [35,36] in DIS and fmnr [37,38] 
in photoproduction are used. The renormalization and the factorization scales are set to /j, r = 
fif = a/<5 2 + 4m^ in DIS and to [i r — /x/ = ^/pf + 4m;? in photoproduction, respectively. For 
both calculations the charm mass is chosen to be m c = 1.5 GeV. The calculations result in 
predictions for the production of charm quarks. To obtain predictions for a measurement of D* 
meson production hadronization corrections evaluated using the LUND hadronization model 
as implemented in RAPGAP are applied. For the longitudinal fragmentation Peterson functions 
are used with e = 0.035 as suggested for NLO predictions by [39]. For the calculation of these 
corrections parton showers are included to simulate the higher order effects of QCD in the event 
generation of the LO Monte Carlo program. To estimate the uncertainty of the NLO calculations 
the renormalization and the factorization scales are simultaneously varied by factors of 1/2 and 
2, the charm mass is varied by ±0.2 GeV and the Peterson fragmentation parameter e is varied 
by ±0.025. The uncertainties originating from all these variations are added in quadrature. 
They result in a combined uncertainty on the theoretical integrated D* meson cross section of 
« 25% in DIS and ss 22% in photoproduction. 

7.2 Two Gluon Exchange Models 

The measured D* cross sections are compared with QCD calculations based on the perturbative 
two gluon approach of 'BJKLW [8] using the k t unintegrated gluon density J2003 set2 evolved 
by the CCFM [40] evolution equations obtained from fits [9] to the inclusive DIS cross section. 
These calculations are applicable only in the region of small xjp (xjp < 0.01), where contri- 
butions from secondary Reggeon exchanges can be neglected. To ensure that the perturbative 
calculations are applicable a cut on the transverse momentum of the gluon of pf > 2.0 GeV for 
the process 7*p — * ccgp is applied. 
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7.3 The MRW Model 



The measurements of the diffractive charm cross section in DIS are also compared with the ap- 
proach of 'MRW [10] which can be considered to be a hybrid of the two approaches described 
in sections ET] and |7]2j The parameters of the input DPDFs were determined from a fit to the 
HI inclusive diffractive data [10]. At low /3, charm is produced via a 'resolved Pomeron' mech- 
anism by BGF-type processes calculated up to order a 2 , as in the approach of section 177X1 At 
high (3, the perturbative two-gluon state participates directly in the hard interaction via 'photon- 
Pomeron' fusion. This 'direct Pomeron' contribution is similar to the ^*p — > ccp contribution 
of the BJKLW model and depends on the square of the gluon distribution of the proton. 



8 Results 

In DIS the integrated cross section of diffractive D*^ production in the kinematic range of 2 < 

Q 2 < 100 GeV 2 , 0.05 < y < 0.7, x F < 0.04, M Y < 1.6 GeV, |t| < 1 GeV 2 , p t (D*) > 2 GeV 
and \i](D*)\ < 1.5 is measured to be 

a(ep -> e/J* ± X'F) D i S = 234 ± 29(stat.) ± 34(syst.) pb, (8) 

which is in good agreement with the measurement in the same kinematic range in the previous 
analysis by HI [11]. 

In photoproduction the integrated ep cross section of diffractive production in the kine- 
matic range of Q 2 < 0.01 GeV 2 , 0.3 < y < 0.65, x P < 0.04, My < 1.6 GeV, |*| < 1 GeV 2 , 
Pt(D*) > 2 GeV and \r](D*)\ < 1.5 is measured to be 

a(ep -> eD^X'Y)^ = 265 ± 50(stat.) ± 41(syst.) pb. (9) 

A comparison of the measured integrated cross sections in DIS and photoproduction with the 
predictions of the NLO calculations for the two sets of HI 2006 DPDFs (Fit A and Fit B) [3] 
is given in table [T] A good agreement between the data cross sections and the NLO QCD 
calculations is observed. 

The D* meson cross section in DIS is also measured differentially as a function of the 
D* kinematic variables Pt{D*) and r)(D*), the DIS kinematic variables y and Q 2 , and the 
diffractive variables Xp, (3 and Zp s . They are listed in table |3] and shown in Figures |3] and 
The data are compared in the figures with the predictions of the NLO QCD calculations. For 
the cross sections as a function of the D* and DIS kinematic variables the predictions for the 
two sets of DPDFS are similar with both providing a good description of the data. For the 
comparison with the diffractive kinematic quantities the differences in the predictions for the 
two DPDFs are larger, with z°p showing the largest sensitivity, where the steeper fall-off of 
the gluon density in Fit B is reproduced. However, within the present experimental errors and 
theoretical uncertainties these differences cannot be resolved. The good description of the NLO 
QCD calculations for all of the D*, DIS and diffractive kinematic distributions supports the 
assumption of QCD factorization, in particular, the compatibility of the gluon density obtained 
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from scaling violations in the inclusive diffractive cross section with that required to describe 
the D* data. 



In photoproduction the D* cross section is shown differentially as a function of the D* 
kinematic variables pt(D*) and r)(D*) and the kinematic variable y in Figure|7]and as a function 
of the diffractive kinematic variables xjp and z°p in Figure [5J The values are given in table |U 
The data are well described by the theoretical predictions within the larger experimental errors 
for photoproduction. As in DIS the largest sensitivity to the different parameterizations of the 
gluon is evident in the Zp s distribution. The shapes of the z°p distribution for the predictions 
in DIS and 7p are compatible which is due to the fact that both kinematic regimes probe the 
diffractive gluon density at a similar scale. 

The good agreement of the NLO QCD predictions with the measured cross sections ob- 
served in DIS and photoproduction, both in shape and normalization, supports the assumption 
that QCD factorization is applicable in both kinematic regimes. A quantity, which is less sen- 
sitive to the input of diffractive parton density functions and theoretical uncertainties is defined 
by 

R DIS = (^mcas/^theo)^ ( 10 ) 

where a mcas and <j theo denote the measured and the predicted integrated cross section for D* 
production. To reduce theoretical uncertainties due to extrapolations from different regions 
in y the cross section in DIS is further restricted to the range of 0.3 < y < 0.65 as for the 
photoproduction measurement. The DIS cross section in this range is shown in table[T] The ratio 
-^ms ls found to be 1.15±0.40(stat.) ±0.09(syst.), with the systematic uncertainty originating 
from the model uncertainty on the (3 distribution in DIS, the fragmentation uncertainties and the 
uncertainties on the Reggeon contribution. The theoretical uncertainty on -R^ s is ±7%. The 
measurement of -R^is snows no evidence for a suppression of the photoproduction component 
although the statistical error of the measurement is large. 

In Figure |5] an additional comparison of both the NLO QCD calculations and of the predic- 
tion from the perturbative two gluon calculation of BJKLW [8] with differential cross sections 
in the range of validity of the two gluon model (xp < 0.01) are shown. The cross sections are 
given in table |5] Within the uncertainties a good agreement between the data and both the NLO 
QCD calculation and the model of BJLKW is observed. For the two gluon calculation in this 
kinematic range the 7*p — > ccgp contribution is seen to dominate with the YP ~^ C ^P process 
contributing only at high values of z°p s . Varying the p t cut-off for the gluon in the 7*p — > ccgp 
process by ±0.5 GeV leads to a variation of the cross section of ~ 25% and is also compatible 
with the data. 

The measurements of the diffractive charm DIS cross sections in Q 2 , y and M\ obtained 
from the displaced track method are converted to measurements in xp, (3 and Q 2 using 

dVg d 3 cr<§ sy 2 



dxjpd(3dQ 2 &Q 2 dydM 2 x (3 
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The diffractive charm reduced cross section is defined as 



5${x P , 0, Q 2 ) - dx ^^ D dQ 2 27Ta 2 {1 %_ y) 2y ( 12 ) 

where a is the fine structure constant. The reduced cross section is approximately equal to the 
charm contribution F®^ cc to the diffractive structure function . The difference is due to 
the contribution from the longitudinal diffractive charm cross section, which is expected to be 
small for the data points presented in this paper. 

The measurements of xjpd c £ obtained from the displaced track method are listed in table 
and shown in Figure [TO] as a function of (5 for fixed values of Q 2 and xjp. In the figure, the 
displaced track method data point measured at x P = 0.01 is interpolated to x P = 0.018 using 
a parameterization of from the NLO QCD fit. The measured points of xjpa^ are compared 
with the results extracted from the D* meson analysis. For this purpose the D* cross section 
is measured in the same Q 2 , y and M\ ranges as for the displaced track method. The results 
are given in table |7J These measurements in the visible D* kinematic range pt(D*) > 2 GeV 
and \rj(D*)\ < 1.5 are extrapolated with the NLO calculation program HVQDIS to the full D* 
kinematic phasespace in order to extract the diffractive open charm cross section. The extrap- 
olation factors are found to be ~ 2.5. The NLO calculation program is also used to evaluate 
the bin center corrections, which are made to the same central values as in the displaced track 
analysis. The HI data are also compared with D* measurements from the ZEUS collabora- 
tion [12] which are interpolated to the same kinematic range as the HI measurement using the 
NLO QCD fit and corrected with a factor of 1.23 to account for the difference in the measured 
range from My = m p to My < 1.6 GeV [24]. The measurements for x pd^ from the displaced 
track analysis and the D* extraction methods from both HI and ZEUS are in good agreement. 
A comparison with the predictions of the NLO DPDFs shows a good description of the data. 

In table [Hand Figure [TT] the measurements are also presented in the form of the fractional 
contribution of charm to the total diffractive ep cross section f^. In the given kinematic range 
the value of /ff is « 20% on average, which is comparable to the charm fraction in the inclu- 
sive cross section at low values of Bjorken x for similar values of Q 2 [14]. The NLO QCD 
predictions shown in Figure [TT] are found to describe the data well. 

In Figures [T2]and[T2]the x F a c ^ and f£ data are compared with the predictions of the MRW 
model [10]. In the kinematic range of the measurements the 'resolved Pomeron' contribution, 
where charm is generated via BGF, is seen to dominate in the model at low j3, while the 'direct 
Pomeron' process, where charm is generated via 'photon-Pomeron' fusion is significant at high 
values of (3. A good description of the data is observed supporting the validity of the DPDFs 
extracted in this model. 



9 Conclusions 



Measurements are presented of the diffractive charm cross section using two independent meth- 
ods of charm reconstruction. In the first method charm quarks are tagged using D* mesons. In 
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the second method tracks, with a significant displacement from the primary vertex, are recon- 
structed using the CST of HI. These displaced tracks arise due to the long lifetime of charmed 
hadrons. 

The diffractive D* cross section is measured in DIS and photoproduction. The integrated 
cross section in DIS is in good agreement with a former measurement of HI, which was ob- 
tained from an independent dataset with less than half the luminosity of the present measure- 
ment. This is the first cross section measurement of diffractive open charm photoproduction 
at HERA. A comparison with QCD calculations in NLO based on DPDFs obtained from 
inclusive diffractive scattering at HI is in good agreement with the measurement in both kine- 
matic regimes. No evidence is observed for a suppression in photoproduction. In the region of 
xp < 0.01 the DIS D* data are found to be also well described by a model based on perturbative 
two gluon exchange and k t - factorization. 

The displaced track measurements are made at Q 2 = 35 GeV 2 for 3 different values of xjp 
and (3. In this kinematic range the charm contribution to the inclusive diffractive cross section 
is found to be « 20% on average which is compatible with the charm fraction in inclusive DIS 
found at low values of Bjorken x for similar values of Q 2 . The cross sections are found to be 
in good agreement with the measurements extrapolated from the D* cross section results and 
to be well described by the predictions of NLO QCD. At low xjp, the data are found to be also 
well described by a hybrid model based on two gluon exchange and diffractive parton densities. 
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HI 99-00 


N(D*) 


Cross Section [ pb ] 








Data 


HI 2006 DPDF 
Fit A Fit B 


DIS 


0.05 < y < 0.7 
0.3 < y <0.65 


122 ± 15 
34 ±8 


234 ± 29 (stat.) ± 34 (syst.) 
55 ± 16 (stat.) ± 9 (syst.) 


287±fJ 
86±?° 


272±^f 
84±?° 


IP 


0.3 < y <0.65 


70 ± 13 


265 ±50 (stat.) ±41 (syst.) 


360±9[] 


359±f 5 



Table 1 : Measured cross sections and NLO QCD predictions for diffractive D* meson produc- 
tion in the visible ranges of DIS and photoproduction (jp). The uncertainty on the NLO QCD 
predictions is given by the variation of the mass, the scale and the fragmentation parameters as 
described in the text. 
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Uncertainty (%) 












(ip) 




iracK ^uio j 


Trigger efficiency 


5 


3 


1 


Scat. e + energy/angle (1% © 1 mrad (DIS), 1.5% ( r yp)) 


2 


2 


2 


Track reconstruction efficiency 


6 


6 


2 


Signai extraction method (D ) 


6 


6 




Ketfections (U ) 


1.5 


1.5 




o resofution (25 /iin © 2UU /im) 






2 


Hadronic energy scale (4%) 


1 


1 


3 


QCD model (reweights in xp, (3, Q 2 , y) 


1 


5 


12 — 18 


Proton diss, model (reweights in \t\, M Y , fraction) 


A 

4 


r 



r 



iNoise in rwiD anu r luu 








Tassins efficiency of FMD (10%) 


1 


1 


1 


Tagging efficiency of PRT (lJo%) 


7 


9 


9 


Non-diffractive background (100%) 


3 


1 


1 


Reggeon contribution (100%) 


4 


1 


1-9 


Fragmentation of c quarks 


4 


1 


7 


Branching fractions / lifetimes / track multiplicities 


2.5 


2.5 


3 


Asymmetry of 5 for light quarks (±50%) 






4-16 


Beauty fraction (tf ° %) 






0-40 


Quark axis (2 °/5 °) 






3 


Luminosity 


1.5 


1.5 


1.5 


QED correction 




2 


2 


Bin center correction 






8-10 


Total 


15 


15 


26-47 



Table 2: Systematic uncertainties for the measurement of diffr active open charm production for 
the inclusive cross section in the visible range for the reconstruction of mesons in DIS and 
photoproduction and in the differential bins for the displaced track method in DIS. 
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DIS D* ± meson cross section as a function of Pt(D*) 


Range (GeV) 


Bin Center (GeV) 


da/dp t (D*) (pb/GeV) 


Sstat 


7c) 


Ssyst (%) 


2.0 


2.5 


2.20 


169 


25 




16 




2.5 


3.0 


2.75 


114 


25 




15 




3.0 


3.6 


3.35 


76 


24 




15 




3.6 


10.0 


5.45 


8 


23 




15 




DIS meson cross section as a function of rj(D*) 


Range 


Bin Center 


da/dr](D*) (pb) 


Sstat (%) 


Ssyst ( 


Jo) 


-1.5 - 


-0.75 


-1.17 


92 


22 




15 




-0.75 - 





-0.33 


101 


21 




15 







0.75 


0.42 


81 


25 




15 




0.75 - 


1.5 


1.12 


40 


37 




21 






DIS D* ± meson cross section as a function of y 


Range 


Bin Center 


da/dy (pb) 


Sstat 


7c) 


Ssyst (' 


7c) 


0.05 - 


0.15 


0.09 


975 


18 




15 




0.15 - 


0.30 


0.22 


423 


25 




15 




0.30 - 


0.45 


0.38 


198 


37 




16 




0.45 - 


0.70 


0.55 


141 


37 




18 




DIS D* ± meson cross section as a function of Q 2 


Range (GeV 2 ) 


Bin Center (GeV 2 ) 


da/dQ 2 (pb/GeV 2 ) 


Sstat 


7c) 


Ssyst ( 


7c) 


2.0 


5.0 


4.0 


17 


27 




17 




5.0 


15.0 


9.5 


7.6 


21 




15 




15.0 - 


35.0 


23.5 


3.6 


21 




14 




35.0 - 


100.0 


60.5 


0.6 


31 




14 




DIS D*^ meson cross section as a function of \og(x]p) 


Range 


Bin Center 


da / dXogixp) (pb) 


Sstat 


7c) 


Ssyst 


7c) 


-3.0 - 


-2.6 


-2.79 


36 


39 




30 




-2.6 - 


-2.2 


-2.39 


118 


22 




21 




-2.2 - 


-1.8 


-2.01 


138 


25 




15 




-1.8 - 


-1.4 


-1.55 


275 


21 




17 




DIS D* ± meson cross section as a function of z°p 


Range 


Bin Center 


da/dzp (pb) 


Sstat (' 


7c) 


Ssyst ( 


7c) 





0.15 


0.07 


312 


44 




19 




0.15 - 


0.45 


0.29 


325 


19 




15 




0.45 - 


1 


0.69 


99 


17 




30 




DIS meson cross section as a function of log 


(P) 


Range 


Bin Center 


da/d\og{(5) (pb) 


Sstat 


7c) 


Ssyst (' 


7c) 


-2.5 - 


-1.8 


-2.12 


55 


40 




19 




-1.8 - 


-1.2 


-1.57 


120 


23 




14 




-1.2 - 


-0.6 


-0.88 


123 


18 




16 




-0.6 - 





-0.28 


50 


27 




21 





Table 3: Differential cross sections for diffractive D* ± meson production in DIS, as a function 
ofp t (D*), r](D*), y, Q 2 , xjp, z°p and (3, given in the range of2<Q 2 < 100 GeV 2 , 0.05 < 
y < 0.7, xp < 0.04, My < 1.6 GeV, \t\ < 1 GeV 2 , p t (D*) > 2 GeV and \r](D*)\ < 1.5. 
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IP 


Z?** meson cross section as a function of pt(D*) 






Range (GeV) 


Bin Center (GeV) 


da/d Pt (D*) (pb/GeV) 


Sstat (%) 


Ssyst (' 


ft) 


2.0 


2.6 


2.25 


160 


39 


16 




2.6 


3.2 


2.95 


172 


26 


14 




3.2 


10.0 


4.95 


10 


29 


17 






IP 


D* ± meson cross section as a function of r)(D*) 






Range 


Bin Center 


da/dr](D*) (pb) 


Sstat (%) 


Ssyst 


ft) 


-1.5 - 


-0.65 


-1.05 


112 


28 


15 




-0.65 - 


0.20 


-0.28 


169 


25 


15 




0.20 - 


1.50 


0.82 


27 


67 


23 




7P meson cross section as a function of y 


Range 


Bin Center 


da/dy (pb) 


Sstat (%) 


Ssyst 


ft) 


0.30 - 


0.40 


0.35 


1010 


34 


16 




0.40 - 


0.50 


0.45 


785 


29 


16 




0.50 - 


0.65 


0.57 


555 


36 


20 




7P meson cross section as a function of log(a;jp) 


Range 


Bin Center 


da/d log (x jp) (pb) 


Sstat (%) 


Ssyst (' 


ft) 


-3.0 - 


-2.2 


-2.59 


77 


30 


19 




-2.2 


-1.8 


-2.01 


266 


27 


15 




-1.8 - 


-1.4 


-1.61 


214 


42 


17 




7P D* ± meson cross section as a function of 


Range 


Bin Center 


da/dzp (pb) 


<Wt (%) 


Ssyst ( 


ft) 





0.15 


0.06 


680 


36 


16 




0.15 - 


0.45 


0.28 


400 


26 


15 




0.45 - 


1 


0.70 


51 


37 


48 





Table 4: Differential cross sections for diffractive D* ± meson production in jp, as a function 
ofp t (D*), rj(D*), y, x F and z°p, given in the range of Q 2 < 0.01 GeV 2 , 0.3 < y < 0.65, 
xjp < 0.04, My < 1.6 GeV, |f| < 1 GeV 2 , p t (D*) > 2 GeV and \r](D*)\ < 1.5. 
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DIS meson cross section as a function of Pt(D*) 


Range (GeV) 


Bin Center (GeV) 


da/d Pt (D*) (pb/GeV) 


Sstat (%) 


Ssyst (%) 


2.0 - 


3.0 


2.45 


58 


23 


18 


3.0 - 


10.0 


4.95 


5 


22 


18 


D 


S meson cross section as a function of r\ 


(D*) 


Range 


Bin Center 


da/dr](D*) (pb) 


Sstat (%) 




-1.5 - 


0. 


-0.65 


48 


18 


17 


0. 


1.5 


-0.33 


16 


34 


23 


DIS D*^ meson cross section as a function of \og(x]p) 


Range 


Bin Center 


da/d\og(x]p) (pb) 


Sstat (%) 


5syst (%) 


-3.0 - 


-2.5 


-2.71 


43 


32 


30 


-2.5 - 


-2.0 


-2.25 


143 


19 


17 


DIS D*^ meson cross section as a function of z°p 


Range 


Bin Center 


da/dzp (pb) 


Sstat (%) 


Sgyst (%) 





0.5 


0.27 


74 


31 


16 


0.5 - 


1. 


0.69 


91 


19 


42 



Table 5: Differential cross sections for diffractive -D meson production in DIS, in the same 
kinematic region as that given in table\3\but further restricted to xjp < 0.01. 



Reduced Cross Section a^(xp, (3, Q 2 ) 


Displaced track 


D* 


Q 2 (GeV 2 ) 




P 


~cc 
D 


Sstat (%) 


5sys (%) 


~cc 
a D 


35 


0.004 


0.25 


1.50 


25 


+27 
-27 


1.33 


35 
35 


0.010 
0.018 


0.10 
0.04 


0.63 
0.62 


23 
18 


+26 
-29 
+29 
-47 


1.20 
0.62 



Table 6: The reduced cross section a^(xp, (3, Q 2 ) obtained from the displaced track method. 
The last column shows the results obtained by extrapolating the D* cross sections in table\7\ 
using the HI NLO QCD fit. 



DIS D* ± meson cross section as a function of Mx 


Range (GeV) 


da/dMx (pb/GeV) 


Sstat (%) 


Ssys (%) 


6-12 


2.5 


45 


20 


12 - 20 


5.0 


26 


15 


20 - 99 


0.39 


42 


17 



Table 7: The differential cross section for diffractive D* production in DIS as a function ofMx 
measured in the range 15 < Q 2 < 100 GeV 2 , 0.07 < y < 0.7, x F < 0.04, M Y < 1.6 GeV, 
\t\ < 1 GeV 2 , p t (D*) > 2.0 GeV and \rj(D*)\ < 1.5. 
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Fractional charm contribution fff 


Displaced track 


D* 


Q 2 (GeV 2 ) 


Xjp 




fee 

Id 


Sstat (%) 


S SyS t (%) 


fee 

Id 


35 


0.004 


0.25 


0.184 


25 


+25 
-25 


0.162 


35 


0.010 


0.10 


0.193 


23 


+23 
-27 


0.367 


35 


0.018 


0.04 


0.278 


18 


+27 
-46 


0.278 



Table 8: The fractional charm contribution to the diffr active cross section f£ obtained from 
the displaced track method. The last column shows the results obtained by extrapolating the 
D* cross sections in table using the HI NLO QCD fit and dividing by the measured total 
diffractive cross section. 
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Figure 2: The AM distribution for each track combination that passes the selections described 
in sections^\and \6.1\ for ( a ) DIS and (b) photoproduction. The parameterization used to obtain 
the number of reconstructed D* mesons shown in the plot is described in the text. 



29 




Figure 3: The significance 5 /a (5) distribution (a) of the highest absolute significance track 
(Si) and (b) of the track with the second highest absolute significance (S2). Included in the 
figure is the expectation from the Monte Carlo simulation program RAPGAP for light, charm 
and beauty quarks. The contributions from the various quark flavors are shown after applying 
the scale factors obtained from the fit to the subtracted significance distributions of the data 
shown in Figure^ 



— 10 3 




Figure 4: The subtracted significance distributions of (a) Si and (b) S 2 . Included in the figure 
is the result from the fit of the Monte Carlo distributions of the various quark flavors to the data. 
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HI Diffractive D* in DIS 
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Figure 5: Differential cross sections for diffractive D* meson production in DIS as a function 
of (a) p t (D*), (b) rj(D*), (c) the inelasticity y and (d) the photon virtuality Q 2 . The inner error 
bars of the data points represent the statistical uncertainties of the measurement only, while 
the outer error bars show the statistical and systematic uncertainties added in quadrature. The 
data are compared with a pQCD calculation in NLO using two alternative sets of diffractive 
parton density functions (Fit A and Fit B) extracted by HI [3]. 
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Figure 6: Differential cross sections for diffractive D* meson production in DIS as a function 
of (a) xjp, (b) z°p and (c) (3. The inner error bars of the data points represent the statistical 
uncertainties of the measurement only, while the outer error bars show the statistical and sys- 
tematic uncertainties added in quadrature. The data are compared with a pQCD calculation in 
NLO using two alternative sets of diffractive parton density functions (Fit A and Fit B) extracted 
by HI [3]. 
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Figure 7: Differential cross sections for diffractive D* meson production in photoproduction 
as a function of (a) Pt(D*) , (b) T)(D*) and (c) the inelasticity y. The inner error bars of 
the data points represent the statistical uncertainties of the measurement only, while the outer 
error bars show the statistical and systematic uncertainties added in quadrature. The data 
are compared with a pQCD calculation in NLO using two alternative sets of diffractive parton 
density functions (Fit A and Fit B) extracted by HI [3]. 
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Figure 8: Differential cross sections for diffractive D* meson production in photoproduction as 
a function of (a) xp and (b) z^ s . The inner error bars of the data points represent the statistical 
uncertainties of the measurement only, while the outer error bars show the statistical and sys- 
tematic uncertainties added in quadrature. The data are compared with a pQCD calculation in 
NLO using two alternative sets of diffractive parton density functions (Fit A and Fit B) extracted 
by HI [3]. 
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Figure 9: Differential cross sections for diffractive D* meson production in DIS, in the re- 
stricted kinematic region ofxjp < 0.01, shown as a function of (a) p t (D*), (b) r)(D*), (c) xp 
and (d) z°p s . The inner error bars of the data points represent the statistical uncertainties of 
the measurement only, while the outer error bars show the statistical and systematic uncertain- 
ties added in quadrature. The data are compared with a pQCD calculation in NLO and to a 
prediction from the perturbative two gluon approach of BJKLW [8] with a cut for the gluon 
momentum in the 7*p — > ccgp process of p t > 2.0 GeV. The dashed line indicates the resolved 
7*p — » ccgp contribution only while the solid line shows the sum of the 7*p — ► ccgp and the 
7*p — > ccp contributions. 
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Figure 10: The measured reduced cross section x p shown as a function of (3 for two different 
values ofxjp. The inner error bars of the data points represent the statistical error, while the 
outer error bars represent the statistical and systematic uncertainties added in quadrature. The 
measurements obtained from D* mesons from HI in this paper and from ZEUS [12] are also 
shown. Measurements at the same values of (3 are displaced for visibility. The measurements 
are compared with NLO predictions based on two alternative sets of diffr active parton density 
functions ( Fit A and Fit B ) extracted by HI [3 ]. 
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Figure 11: The contribution of charm quarks to the total diffractive cross section f£ shown as 
a function of (3 for two different values of xp. The inner error bars of the data points represent 
the statistical uncertainties, while the outer error bars represent the statistical and systematic 
uncertainties added in quadrature. The measurements are compared with NLO predictions 
based on two alternative sets of diffractive parton density functions (Fit A and Fit B) extracted 
by HI [3]. 
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Figure 12: The measured reduced cross section x jp a c £ shown as a function of (3 for two different 
values of xp. The inner error bars of the data points represent the statistical error, while the 
outer error bars represent the statistical and systematic uncertainties added in quadrature. The 
measurements obtained from D* mesons from HI in this paper and from ZEUS [12] are also 
shown. The measurements are compared with the model of MRW [10] based on perturbative 
two gluon exchange and DPDFs. 
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Figure 13: The contribution of charm quarks to the total diffr active cross section shown 
as a function of j3 for two different values of x p. The inner error bars of the data points 
represent the statistical uncertainties, while the outer error bars represent the statistical and 
systematic uncertainties added in quadrature. The measurements are compared with the model 
ofMRW [10] based on perturbative two gluon exchange and DPDFs. 
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